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ABSTRACT 

Carbon Monoxide is a commonly used IR/sub-mm tracer of gas in proto- 
planetary disks. We present an analysis of ultraviolet CO emission in HST-COS 
spectra for 12 Classical T Tauri stars. Several ro-vibrational bands of the CO 
A 1 ]! - A^E 4 " (Fourth Positive) electronic transition system are spectrally resolved 
from emission of other atoms and H 2 . The CO A l U v '=14 state is populated by 
absorption of Lya photons, created at the accretion column on the stellar sur- 
face. For targets with strong CO emission, we model the Lya radiation field as 
an input for a simple fluorescence model to estimate CO rotational excitation 
temperatures and column densities. Typical column densities range from Nco — 
10 18 - 10 19 cm" 2 . Our measured excitation temperatures are mostly below Tqo 
= 600 K, cooler than typical M-band CO emission. These temperatures and the 
emission line widths imply that the UV emission originates in a different popula- 
tion of CO than that which is IR-emitting. We also find a significant correlation 
between CO emission and the disk accretion rate M acc and age. Our analysis 
shows that ultraviolet CO emission can be a useful diagnostic of CTTS disk gas. 

Subject headings: protoplanetary disks — stars: individual (RECX 15, DM Tau, 
DE Tau, HD 135344B, DR Tau, DK Tau, DF Tau, LkCa 15, GM Aur, HN Tau, 
RECX 11, V4046 Sgr) 
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Introduction 



Classical T Tauri stars (CTTSs) provide the opportunity to witness the birth of extra- 
solar planetar y systems. Their disks of ga s and dust are the environments in which Jovian 
planets form (jPullemond fc Monnierll2010l ). a process that we are only beginning to under- 
stand. From the ground, carbon monoxide (CO) is the easiest molecule to obse rve in emis- 



sion f rom these disks, and has been studied from near-IR to mm wavelengths ( iSalyk et al. 



20081 : iDutrev et al. 



2008 



Sturm et al.l 120101 ). Such observatio ns sample regions from r in < 



1 AU (INaiita et al.ll2003h to r out > 800 AU flPietu et al.ll2007t ). The second most abundant 
molecule in the galaxy (after H 2 ), CO plays a significant role in the formation of planetary 
systems. As H 2 does not emit strongly in the IR, CO is an important tracer of disk gas. 
CO abundance varies greatly in different environments throughout the M ilky Way. In the 
diffuse ISM, the CO/H 2 ratio ran ges from 10~ 7 to 10~ 6 (IBurgh et al.l 120071 ); in dense clouds, 
i t can be 10~ 4 (ILacv et al.l Il994f) , whil e comets in the solar system exhibit a ratio of ~10 
(ILupu et al.l 120071 : iFeldman et al.ll2009l ). During the process of star and planet formation, 
the CO and H 2 may be s electively removed from or added to the gas p hase in a way that 
changes the CO/H 2 ratio (jVidal-Madjar et al.lll994j : iFrance et al.ll2011ah . By probing CO in 
protoplanetary disks, we can constrain models of planetary system formation. 

Electronic transitions of CO have recently been identified in the H ubble Space Telescop e 
Cosmic Origins Spectrograph (HST-COS) spectra of several CTTSs (IFrance et al.l 1201 lah . 
Several vibrational bands of the A 1 !! - V 1 £ + (Fourth Positive) system were measured in 
emission for V4046 Sgr and RECX 11, and in absorption in HN Tau. Lyman a (Lya) 
photons produced in the accretion funnel flow excite the CO, populating its A 1 !! v' = 14 
state; several of the resulting (14 - 2) through (14 - 12) transitions are detected between 
1280 A and 1720 A. Emission from the (0 - 1) 1600 A and (0 - 3) 1713 A bands also appears, 
potentially due to C IV photoexcitation of CO. The shape and strength of these emission 
bands vary with temperature (Tco) and column density (Nco) such that these parameters 
can be determined. Circumstellar CO can be observed in absorption as well; a collisionally- 
excited distribution of gas at T C o = 500 ± 200 K is able to reproduce several of the (v' - 0) 
absorption bands in the HN Tau spectrum for logio(iVc'o) = 17.3 ± 0.3 cm -2 . France et al. 
(2011, in press) also find CO in absorption in AA Tau, with \og\o(Nco) = 17.5 ± 0.5 cm -2 . 

We extend this work by conducting a survey of COS spectra for 12 CTTSs in greater 
detail. Our aim in this paper is to characterize T C o and N C o m a variety of CTTSs, so as 
to understand correlations between the UV-emitting CO (UV-CO) and the disk properties 
that characterize these systems. The UV-CO can also be compared with CO previously 
detected in these disks at IR and mm wavelengths. We model the Lya radiation field 
assuming absorption from an outflow and interstellar H I, using this as an input for simple 
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fluorescence models. These models enable us to estimate Tco and Nco for targets with high 
signal-to-noise (S/N) CO emission. In Sj2]we detail the COS observations, then discuss the 
CO emission in section §21 In §H we summarize our model for the Lya profiles as well as 
the CO fluorescence model. In §4.21 we discuss the results of this analysis, and our ability 
to constrain the temperatures and column densities of CO in each of the targets. We also 
compare these empirically derived properties with other previously measured parameters 
from the literature. 



2. Observations 



We investigated a sample of 12 CTTSs observed by H ST-COS; the observational pa- 
rameters are listed in Table [I] in order of observation date. These were the first CTTSs 
observed with COS. Spectra from several of these targets have already been presented in 
the literature. Ilnglebv et al.l ( 1201 ll ) studied accretion rates in RECX 11 and RECX 1, and 
Yang et al.l (120111) demonst r ated L ya photon absorption and re-emission by H 2 in V4046 Sgr 



and DF Tau. IFrance et al.l (2011b ) characterized the continuum emission of V4046 Sgr and 
DF Tau, while IFrance et al.l ( l2011af ) identified CO Fourth Positive emission and absorption 
in the spectr a of V4 046 Sgr, RECX 11, and HN Tau. This paper builds upon the work of 
France et al.l (l2011a| ). advancing the models for analyzing the CO emission of a survey of 
CTTSs. 



Each spectrum consists of 2 - 4 orbits with the G130M and G160M gratings (lOsterman et al 



201 ll ). yielding a moderate resolving power of X/5X ~18,000. As there is no field stop in the 
instrument, the entire objects were sampled without appreciable spatial resolution relative 
to the size of the disks. We used CALCOS, the COS calibration pipeline, to reduce the 
data and then created coadded spectra from 1140 A to 1760 A. The data red uction pipeline 



creat es inaccurately small error bars in low count regions of the spectrum ( IFroning et al. 



201 ll ). This gives excessive weight to continuum points in the least squares fits of the emis- 
sion feature. We empirically determined that this was generally the case for data points with 
errors more than three times smaller than the mean la continuum error. Any errors smaller 
than this limit were thus set to the limit value. 

Most targets belong to the rj Chamaeleontis and Taurus-Auriga star formation regions, 
with properties described in Appendix |A] Table |5]lists a variety of target parameters relevant 
to this work, including visual extinction (Ay), accretion rate (M acc ), inclination, and age. 
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3. Lya-excited CO Emission 

CO UV emission appears in 7 of the 12 targets in our survey. We find strong CO (14 - 
v") emission in four targets (V4046 Sgr, DF Tau, RECX 15, and DM Tau), and weaker but 
identifiable emission in RECX 11, HD 135344 B, and DE Tau. The CO emission bands seen 
are listed in Table EJ The (14-1), (14-6), (14-9), and (14 - 11) bands are unidentifiable due 
to either low branching ratios or overlying emission from other molecules and ions. We list 
the branching ratios B(v' - v") for each band, calculated by summing the A values of each J 
transition in a given band, and dividing by the sum of all (v' - v") transitions for a given v '. 
Despite the relatively low A values of the overlying (14 - 0) transitions (1.32 xlO 4 s _1 for J" 
= 0), there is sufficient Lya flux to yield strong CO emission. C IV-pumped emission bands 
are also seen in several targets; we present this evidence in §5.31 The Lya-excited emission 
is the primary subject of our analysis. 

Figures [T] and [2] show the observed Lya profiles and (14 - 3) CO bands for the sample. In 
this paper we use femto-erg flux units (FEFU, where 1 FEFU = lxlO -15 ergs cm -2 s~ x A -1 ) 
as a standard flux unit. The geocoronal Lya emission in the middle of the Lya profile is set 
to zero to avoid confusion with target Lya emission. For future discussion, we refer to the 
Lya emission on the longer wavelength side of the interstellar absorption trough as the "red 
wing", and the emission on the shorter wavelength side the "blue wing". The CO plots in 
these figures also separate the CO bands into a "high- J" and "low- J" side, identified by the 
color of the continuum. The blueward CO emission comes from the lower- J transitions (0 - 
15), which are closely spaced relative to the instrument resolution, forming a "bandhead". 
The redward CO emission is composed of higher J transitions, which are well separated in 
several targets. There is a general correlation between the relative amount of flux in the 
Lya blue and red wing and the relative amount of flux seen in the low- J and high- J CO 
emission. For example, in RECX 15 and DF Tau, weak blue wing Lya emission correlates 
with the lack of low- J emission in their CO bands, whereas in V4046 Sgr we observe strong 
blue and red emission in both the Lya and CO profiles. This will motivate our modeling, as 
discussed in §U 

The (14 - 3) band is mostly free of emission from other molecules or atoms; only the 
weak H2 (4 - 2) P(19) line at 1316.55 A appears in all targets. There is an H2 (3 - 2) P(18) 
line at 1320.22 A, but it appears strongly only in DF Tau. Overall, the low- J bandhead 
and separated high-J lines of the (14 - 3) band are uncontaminated, which is crucial for 
determining Nco an d Tco- F° r this reason, we primarily use the (14 - 3) band for analysis 
in all targets. By branching ratios, the (14 - 4) band is stronger, however strong 01 ] (1355.71 
A) emission lines appear over the high-J portion of the band in most targets. Figure [3] shows 
the (14 - 4) band for our sample, with similar continuum denotations as the (14 - 3) band. 
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Only V4046 Sgr, RECX 15, and DM Tau are free from contamination on the high- J side. 
For these three targets we analyze both the (14 - 3) and (14 - 4) bands. 

In V4046 Sgr, DF Tau, DM Tau, and RECX 15, the high-J lines are well separated. 
The high-J lines of DF Tau, RECX 15, and DM Tau appear at the width of the COS line- 
spread-function. In RECX 11 and HD 135344B the CO high-J lines are not distinguishable 
from one another, although whether this is physical (due to velocity broadening of close-in 
CO) or merely due to the lower S/N in the band is not clear. The Lya profile of DE Tau 
is highly absorbed, yet the (14 - 3) features are evident to the eye. HN Tau, DR Tau, DK 
Tau, LkCa 15, and GM Aur show little CO emission above the continuum. Of those, GM 
Aur is the only target to have a Lya profile that is mostly unabsorbed by circumstellar and 
interstellar H I. 



4. Method and Modeling Results 
4.1. Lya and CO Fluorescence Model 



France et al.l (l2011al ) utilized a simple fluoresence model to calculate the temperature 
T C o an d column density N C o for two CTTSs. However, they used the observed Lya emission 
to excite the CO. In our full sample, the targets demonstrate the need for a more detailed 
estimation of the incident Lya radiation; using the observed Lya profiles as the incident 
radiation field cannot reproduce the observed CO band shape. For example, with DF Tau, 
more red wing Lya emission than is seen is necessary to yield the strong CO (14-3) high-J 
line emission observed, relative to its weak low- J emission. We therefore introduce a straight- 
forward model of the physical situation, shown in Figure HJ which we use to re-construct the 
incident radiation. Lya emission is created at the accretion funnel flow near the proto-stellar 
surface, as well as in the photosphere itself. The plot in Figure H] shows such a model Lya 
profile for V4046 Sgr in red. Along our line of sight the Lya is then absorbed by H I in an 
outflow (magenta), followed by interstellar gas (green). Our method consists of two steps: 
first, we create a series of model Lya profiles that fit the observed profile after absorption 
by outflow and ISM H I. Second, we use each outflow-absorbed model Lya emission in a CO 
fluorescence model to match the observed CO emission. 

The observed Lya emission is modeled as: 

1 obs ^ Lya ^Lya^ SM) (1 ) 

where Ciya is the continuum intensity near Lya and the source Lya profile Ii ya is approx- 
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imated as the sum of two Gaussian emission components, parameterized by peaks I\ and 
I2, as well as widths o\ and a 2 . The first Gaussian is at rest relative to the star, while the 
second Gaussian component is centered at a velocity v 2 - The outflow absorber r out is charac- 
terized by a Voigt profile with N out at a velocity v out , and is necessary to match blue-shifted 
absorption seen in most of the targets. The ISM absorber tjsm is determined by a Voigt 
profile with N out at the local standard of rest (LSR). Ay values from Tabl e [2] and Ry 



4.0 we re used to de-redden the data, according to the extinction curve from ICardelli et al. 



(119891 ) . Seven parameters determine the two emission (stellar and accretion shock) and two 
absorption (interstellar and outflow) components. As a result, many parameter combina- 
tions exist that in the end fit the observed Lya profile. However, only some will adequately 
produce the observed CO emission, once run through the fluorescence model. 



Our CO fluorescence model uses a method similar to that implemented by iFrance et al. 



(l2011aT ) . The g round rotational ley els are populated thermally with Tco, using Dunham 



coefficients from iGeorge et al.l (119941 ) to calculate ground state energy values. The Lya flux 
absorbed by the (14 - 0) band is re-emitted according to the branching ratios for the (14 
- v") transitions. However, instead of using the observed Lya flux, we excite the CO with 
the outflow-absorbed profile. For each Lya profile produced in the first step, we produce a 
model CO emission band to fit to the data. The variable parameters in the CO model are the 
gas rotational temperature (Tco), column density (Nco), and the velocity broadening of the 
emitting gas (Avb roa d). This Avb roa d value is the full width at half maximum (FWHM) of the 
Gaussian emission line, assuming the broadening is primarily due to orbital motion. For a 
continuum, we use the average continuum level from nearby areas of the spectrum. Assuming 
that the electronic ground state population can be described by a thermalize d population, 



the res ulting CO emission is primarily dependent on Nco an d Tco- As in France et al. 



( l2011af ). we make no assumptions about the spatial distribution of the CO around the star. 
Due to the low transition probabilities of the (14 - 0) band (~ 10 4 s" 1 ), optical depth effects 
are minimal; therefore, Nqo an d the filling fraction of the CO are degenerate. 

There are limitations to our assumption of a single CO rotational temperature, as mul- 
tiple temperatures are possible throughout the disk. The gas disk temperature varies based 
on height in the disk and distance from the star. As the entire disk is sampled in the COS 
aperture, CO may be observed at multiple temperatures. However, a single temperature 
model does enable us to determine general characteristics of the CO in a large sample with 
minimal parameters. Future models of individual targets should take into account the effects 
of density and temperature on the radiative transfer of the Lya radiation and fluorescent 
emission. 
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4.2. Results 

In half of the positive detection targets, the uncertainties on Nco are less than 0.5 dex 
and the uncertainties on Tqo are less than ~200 K. High S/N in both the Lya and CO 
profile are required for this. In most of the positive detections, values for logio(iV co ) are 
between 18 and 19 and temperatures lie between 100 and 600 K. 



4.2.1. Lya Model Fits 

The top panel of Figure [5] shows representative Lya model results for DF Tau. Several 
model Lya profiles are plotted over the observed Lya profile. The red, blue, and green 
lines correspond to the maximum, average, and minimum Lya emission profiles after being 
absorbed by the outflow H I of each respective model. These are displayed to convey the 
variety of incident radiation profiles that, when used in the CO fluorescence model, produce 
CO emission that matches the data. The magenta line shows a typical outflow + ISM- 
absorbed profile that fits the Lya profile. In the case of DF Tau the entire blue wing is 
absorbed, which we attribute to a non-gaussian absorber velocity profile; therefore we only 
fit the red wing. The blue wing is generally fitted between 1210 A and 1214 A, while the 
red wing is fitted from 1217 A to 1220 A. We did not use two Lya emission components for 
DF Tau and RECX 11, as a second component did not significantly improve the fit. Strong 
H9 absorption ex ists in the Lya wings of V4046 Sgr, DF Tau, RECX 15, and DM Tau 



( lYang et al.ll201lf ); a spline interpolation of the profile was used to minimize their influence 
on the fits. The H 2 absorption may be present in the Lya profile incident on the fluoresced 
CO, but the absorption has a negligible effect on the resulting CO fluorescence. 

We searched a grid of the parameters to find the best-fitting combinations. A "good" 
fit was determined by minimizing \ 2 - Fits were excluded with x 2 greater than 95.4 % (2a) 
from the peak of the \ 2 probability distribution. In Tables H] and we list parameter values 
of the best-fitting Lya profiles. Since our model Lya profiles must not only fit the ob- 
served Lya emission but also produce the appropriate CO emission through our fluorescence 
model, we only include models that create CO emission that fits the (14 - 3) band. The 
Lya parameters listed are determined by taking the average and standard deviation of each 
parameter for all models. The broad range of values for the emission components reflects 
the variety of emission profiles that in combination with different absorbing profiles match 
the data. Continuum values are determined by averaging emission-free ranges of the sur- 
rounding spectrum. The "outflow" absorber velocites for DM Tau are positive, possibly due 
to predominantly inflowing material. Overall, o ur combined values for Nt$m and N nuf are 



similar to those measured previously in CTTSs (IWalter et al.ll2003l ; iHerczeg et al.l 120061 ) . 
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4-2.2. CO Model Fits 

The middle panel of Figure [5] shows a typical CO model fit to the (14 - 3) band, along 
with its value for Nco and Tco- However, this Nco and Tco combination is the best fit for 
only one of our model Lya profiles. The bottom panel shows the distribution of best fit Tqo 
vs. Nco f° r all the best fit Lya profiles, to demonstrate the range of possibilities enabled 
by our model. Each point corresponds to a unique set of Lya profile, Nco, and Tco that in 
combination adequately reproduce the observed Lya and CO emission. The error bars for 
each fit are not included, as they are smaller than the overall distribution of points. Lya and 
CO must be significantly detected in order to sufficiently constrain Nco and Tco', even for 
the targets with the highest quality Lya and CO profiles (V4046 Sgr, DF Tau, and RECX 
15), the uncertainty in the Lya incident radiation causes a spread in the best-fitting Nco 
and T C o- 

For fitting the (14 - 3) band we choose the regions between 1315 A and 1322 A carefully 
so as to avoid the 1316.55 and 1320.22 A H 2 lines. We include the bandhead (J" < 5), as it 
constrains Nco for low Tco- Also, the high J lines 1317 to 1320 A must be included, as they 
constrain Tco for a given value of Nco- The range between 1316.5 and 1317.5 A is generally 
excluded; interstellar H I absorption removes any knowledge of the Lya flux which results 
in CO emission in this region. This corresponds to the J'=6-9, 9-12, and 12-15 transitions 
for the P, Q, and R branches, respectively. For the (14 - 4) band in V4046 Sgr, RECX 15, 
and DM Tau, we exclude flux between 1353.5 A and 1354.5 for the same reason. Due to 
the computation time necessary to perform a grid search with our fluorescence model, we 
utilized a parameter optimization program starting from an initial coarse parameter grid. 

Table [6] lists the CO fit parameters for all targets. Individual confidence intervals for 
each fit are smaller than the overall distribution of good fits. For this reason, the CO 
parameters we list are determined from the minimum and maximum parameter values, with 
the average of the two chosen as the mean value. Information for the (14-4) fits is also 
included for the three targets free of O I ] emission (V4046 Sgr, RECX 15, and DM Tau). 
For targets with both (14 - 3) and (14-4) bands fits, the values of T C o, N C o, and At> broad 
agree within error bars. The required vco to fit each band differs by <10 km s -1 , which is 
below the accuracy of the COS wavelength solution. The lack of knowledge about the core 
of the Lya profile is chiefly responsible for discrepancies between the model flux and that 
of the observed spectrum. Flux in the 1316 - 1317 A region is typically under-predicted, 
in part due to the contaminating H2 line, but also due to a more complicated absorbing 
outflow velocity distribution than we can reasonably model at present. Future work will 
include more complicated velocity profiles. 

In DF Tau, RECX 15, and DM Tau, the high-J lines appear close to unresolved, which 
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help to strongly constrain Avbroad- For DM Tau, HD 135344B, and DE Tau, the Lya and/or 
CO S/N is low enough to make it difficult to constrain parameters as well as in the other 
targets. In RECX 11 the CO high- J lines are not resolved, although whether this is physical 
(due to velocity broadening of close-in CO) or merely due to the lower S/N in the band 
is not clear. As a result, Avbroad has a wide range of values. DE Tau has a wide range 
of Tco an d Nco, primarily due to its poorly-constrained Lya. In GM Aur, HN Tau, DR 
Tau, DK Tau, and LkCa 15, the highly-absorbed Lya and/or weak CO emission makes it 
impossible to constrain N C o and T co better than several dex. Upper limits on N C o cannot 
even be determined without better knowledge of their Lya profiles. This will be possible 
in future work by partially reconstructing the Lya profiles with measured H 2 fluxes. The 
models presented here assume isothermal CO. Relaxing this assumption would provide more 
flexibility in simultaneously fitting low-J and high-J emission lines. The COS observations 
of some TTS may have enough sensitivity to test disk models that predict a range of gas 
temperatures. 



Discussion 



5.1. Spatial Distribution of UV-CO 



While the degeneracy of Nco and the filling fraction prevents us from definitively deter- 
mining the geometry of the UV-CO around the star, comparing our measured UV-CO prop- 
erties with those of previously measured disk gas still reveals much about the spatial distribu- 
tion of the gas. C O M-band fundamental emission [v'=l, v"=0) has been detected in several 



of th ese targets (jSalyk et al.l 120091 ; iNajita et al.l 120031 ; iBast et al.l l201ll ; iPontoppidan et al. 



20081 ). The gas in which the IR-CO emission originates is typically warm, thermally excited 
at temperatures as high as ~1500 K. Excitation diagrams of larger IR- CO samples can b e 
adequately described by LTE models of a single-temperature slab of gas (jSalyk et al.ll2011bl ). 
Line-width analysis suggests that this emission is produced < 1 AU from the central star. 
This information, combined with our measured UV-CO properties, presents two arguments 
that the UV-CO and IR-CO are spatially separate. 

The UV-CO temperatures suggest separate populations, as our measured T C o values 
are several times lower than typical IR-CO temperatures. The ro-vibrational levels of the 
CO ground electronic state are primarily populated based on the Einstein A coefficients of 
each v and J level, relative to the temperature-dependant collision rates and gas density. 
Models typically assume the IR-CO lies in a "skin" of the gas disk where de nsities are high 



enou gh for the rotational temperature to be equal to the kinetic temperature (IBrittain et al. 
20071 ). If this is also true for the UV-CO, then it must lie in a different region of the disk 
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where the gas kinetic temperature is lower. However, as temperature and density vary based 
on location in the disk it may be possible to find gas populations with thermal low- J levels, 
while the high- J levels are subthermal. In targets with well-defined high- J lines (V4046 Sgr, 
DF Tau, RECX 15, and DM Tau), these lines strongly influence the Tqo and Nco fits- If 
densities are less than the critical densities for these high-J states (n cr ji for J =18 at 400 
K is ~ 10 6 cm -3 ), then the derived temperature could be less than that determined from 
the thermalized lower states. In this case the CO density would have to be lower than the 
thermalized IR-CO, and therefore be in a different region of the disk. 

Our measured Avbroad values enable us to spatially constrain the UV-CO for the targets 
with well-defined high-J lines. Due to the blending of the UV-CO bands we cannot calculate 
the line width at zer o intensity of individ ual J transitions, from which the inner CO radius is 



typically calculated ( iBrittain et al.ll2009l ). However, assuming that the FWHM of the high-J 
lines is primarily due to orbital motion of the gas, we can calculate an average orbital radius 
Rco- 



CO 



2sin{i) 
Avb, 



road 



(2) 



where G is the gravitational constant, M* the stellar mass, and i the disk inclination, pre- 
viously measured in the literature. For the values listed in Appendix |A] we calculate Rco 
for V4046 Sgr, DF Tau, RECX 15, and DM Tau: 2.76 AU, 2.16 AU, >3.04 AU, and >1.56 
AU respectively. These radii are outside typical IR-CO radii { Rco < 1 AU) calcu lated with 
Equation [2] and FWHM values from a survey of CTTSs from ISalyk et al.l (j2011al ) . DF Tau 
is observed in t heir survey, with Rco = 0.15 AU based on their CO FWHM (79 km s" 1 ), 
inclination (85°, I.Iohns-Krull fc Valentil ( boOlh ). and M* (0.27 M , ICullbring et~all jl998h ). 
The Avbroad values of RECX 15 and DM Tau are at the COS resolution, so we can only 
determine lower limits for Rco i n these disks. If thermal broadening, turbulent velocity, 
or disk winds contribute to the line widths, then Rco could be even larger. Regardless, 
the temperature and Avbroad disagreements lead us to the conclusion that the UV-CO and 
IR-CO are indeed separate populations. 

In addition to the spatial differences between the two CO populations, there is no con- 
sistent agreement between the detec tion of IR-CO and UV-CO emission in our sample. DR 
Tau exhibits strong IR-CO emission ( ISalyk et al.ll2008f ). but no UV-CO emission is observed. 
Conversely, DM Tau shows strong UV-CO emission, but no IR-CO emission. DF Tau, HD 
135344 B, and DE Tau show both. Unfortunately, the stronger detection targets V4046 Sgr, 
RECX 15, and RECX 11 have yet to be observed in the M-band. Different morphology 
(UV-CO only, IR-CO only, both) is not a strong argument for separate gas populations. In 
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principle, different morphologies could occur in a single gas population if the Lya pumping 
and IR excitation mechanisms are not physically related. Further study of a larger survey 
will be necessary to draw any solid conclusions on this matter, but the temperature and 
radius differences found in our analysis argue against a single gas population. 

The UV-CO also appears to differ from H 2 previously studied in CTTSs. H 2 is best 
observed in the UV due to its lack of a permanent dipole moment. In the TW Hya disk, 



Herczeg et al.l (120041 ) detected H 2 emission, measuring T H2 = 250015QQ K and \ogN H2 = 
18.5^Qg cm -2 . In addition, France et al. (2011, in press) observed H 2 in absorption in the 
spectrum of AA Tau, measuring T H2 = 2500^ K and \ogN H2 = 17.9±g;| cm" 2 . As with the 
IR-CO, the differences in temperature suggest that the two populations reside in different 
regions of the disk. Similar to the conclusion of France et al. (2011, in press) about the 
AA Tau CO and H 2 , a subthermal CO population could explain the disparity. If the CO 
and H 2 are indeed cospatial, then the similarit y of Nco and Nh o derived from UV spectra 



would seem to violate the CO/H 2 ratio of 10 4 ( jLacy et al.lll994l ). It is possible that during 



the process of star and planet formation, the CO is enhanced and/or H 2 removed f rom th e 



gas phase in a way that increases the CO/H 2 ratio. As suggested by France et al.l (j2011al ). 
collisions of large grains or planetesimals may re-release gas-phase CO into the disk. In 
addition, H 2 could be selectively dissociated by radiation in the inner disk. However, in 
order to truly compare both molecules the UV H 2 and CO emission must be studied with a 
self-consistent approach. 

Comparing N co and T C o with other disk properties published in the literature does 
not yield any significant correlations. For example, an increase in rotational temperature vs. 
accretion rate might be expected; higher accretion rates create greater FUV continuum and 
Lya emission, which in turn heats the gas to higher temperatures via photo-electric heating. 
The FUV photons could also populate higher J states in a way not possible by collisional 
processes alone. However, as seen in Figure EJ there appears to be no correlation between 
temperature and accretion rate. Over more than two orders of magnitude of accretion rate, 
targets with strong CO and Lya emission well constrain the temperature between 200 K 
and 500 K with no apparent trend. Only in DM Tau and DE Tau is the temperature less 
constrained, but this is primarily due to their highly absorbed Lya profiles. 



5.2. CO Band Flux and Disk Evolution 

The total CO emission from our sample can be evaluated by integrating the (14 - 3) band 
flux. Referring to Figures [I] and EJ the (14 - 3) bands are summed between the dashed grey 
lines at 1315.2 A and 1320 A. The continuum level C(i4_3) is determined from the surrounding 
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line-free spectrum, and denoted by the horizontal solid black/blue/red line. The blue and 
red lines identify the low and high J state emission regions, which are separately integrated. 
The dashed orange lines above and below mark the 1-a error of the continuum level. The 
integrated flux values are listed in Table UJ For DR Tau and LkCa 15 we list zero flux, 
which actually corresponds to an overall negative flux. Shown in Figure UJ is ratio of Lya 
blue to red flux vs. the ratio of CO blue (low- J) to red (high- J) flux. This demonstrates 
our initial correlation between the relative strength of the Lya photo-exciting radiation and 
the resulting CO emission. 

To test correlations with disk evolutionary properties we use the total CO (14 - v") 
luminosity, calculated from the (14 - 3) band fluxes: 



L{CO) u ^ 



F(14-3) 
5(14-3) 



And 2 



(3) 



where F(14 - 3) is the integrated band flux, B(14 - 3) is the branching ratio listed in Tabled 
and d is the distance to the target. In Figure [BJ L{CO)n-^ shows a strong correlation with 
mass accretion rate. This is not surprising, as the amount of photo-exciting Lya is expected 
to increase with accretion. Combined with the weak opacity of the CO (14 - 0) band, we are 
simply seeing the re-emission of most of the Lya photons capable of being absorbed. 

As a result, M acc cannot be used as a pr oper evolutionary diagnostic, despite th e fact 
that accretion does decrease as a CTTS ages (IHartmann et al.lll998l ; iFang et al.l 120091 ). We 
attempt to correct for this by dividing the CO luminosity by M acc to compare with age. 
Figure [9] shows this relation, with a trend towards more CO emission in more evolved systems. 
One interpretation is that as inner disks clear, the Lya can more readily illuminate the CO. 
V4046 Sgr, HD135344B, DM Tau, and GM Aur are all confirmed transitional objects with 



cleared dust holes i n the center of their disks ( Uensen fc Mathieulll997l ; iRodriguez et al.ll2010 



Brown et al.l l2009r ICalvet et al.l 120051 ) , which would allow the Lya photons to reach more 
CO. CO-IR emission from HD 135344B and GM Aur support t his argument, as the CO- 
IR appears inside the inner dust radii ( Pontoppidan et all 2008 ). However, no UV-CO is 
detected in the transitional disk LkCa 15 A second inte rpretation would be the increase of 
gas-phase CO through grain and planetesimal collisions (IFrance et al.ll2011al ). CO produced 
in the dust disk could be transported to the inner disk surface as the system evolves. 
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5.3. C IV Photoexcitation 

Finally, we measure flux in the 1712 - 1720 A region, assuming that emission in this 
region is due to the (14 - 12) band. Of particular note is how large these fluxes are compared 
to what would be expected from simply applying branching ratios to the (14 - 3) fluxes. 
Figure [TD] shows this graphically; for one of the best-fit (14 - 3) models of RECX 11, the 
corresponding (14 - 12) emission is plotted over its observed 1712 A emission. This trend is 
seen significantly in all targets and also shown in the bottom of Figure ITUl with the dashed line 
representing the expected (14 - 12) integrated flux for a given (14 - 3) integrated flux. If this 
emission is due to CO, some process other than Lya photo-excitation must be responsible. 
It is likely the result of C IV photo-excitation, which is capable of exciting the v'=0 level of 
CO fjFrance et al. 2011a ). The resulting (0 - 3) transition creates a band of emission at 1712 



- 1720 A. Cursory analysis shows that there is sufficient flux in the observed C IV lines to 
produce these levels of emission. 



6. Conclusions 

In this paper, we have characterized UV photo-excited CO emission in a sample of T 
Tauri star spectra from the COS instrument on HST, reaching the following conclusions: 

1. Electronic transitions of CO are observed in approximately half of the T Tauri stars 
surveyed, primarily due to photo-excitation by Lya emission coincident with the A l U - V X S + 
(14 - 0) absorption band. The resulting transitions from the v' = 14 state appear as emission 
bands between 1280 A and 1712 A. The (14 - 3) and (14 - 4) bands are most conducive to 
spectral modeling. 

2. The observed Lya emission is not sufficient to reproduce the observed CO band 
emission; thus, we create a model to simulate the local Lya radiation field incident on the 
CO molecules. With two Gaussian emission components created at the star, followed by 
outflow and ISM absorbers, we are able to recreate the observed Lya profiles and CO band 
emission. 

3. For seven of the twelve targets, we measure CO temperatures between ~100 and 
1500 K, and CO column densities of ~ 10 18 - 10 19 cm -2 . Both strong Lya and CO emission 
are required in order to constrain Nqo to within 0.5 dex; for these 5 targets temperatures 
are below 600 K. The high- J lines in several targets are narrow enough to constrain the CO 
to outside 1 AU. 

4. There are no strong correlations between our calculated Nco and Too values and 
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other previously determined parameters. Our measured CO temperatures are lower than 
previously measured IR-CO and UV-H 2 in CTTSs. The temperature differences, along 
with our calculated emission line widths, imply that the UV-CO and IR-CO are different 
populations. 

5. If the UV-CO and UV-H2 are co-spatial, our Nco values are 3-4 orders of magni- 
tude greater than those expected from previously measured H 2 abundances and dense cloud 
CO/H 2 ratios. Further analysis of the H 2 lines, as well as more detailed geometric modeling, 
is essential to understand the relationship between the UV-CO, IR-CO and H 2 

6. CO emission increases with M acc because of the corresponding increase in the Lya 
radiation field. Upon correcting for M acc , there is still a correlation between CO emission 
and age. We conclude that as CTTSs evolve and disk dust and H I clears, the Lya photons 
are able to reach more CO. 
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NNX08AC14G, NNX08AX05H, and NAS5-98043 to the University of Colorado as well as 
STScI grants to program GO-11616 and GO-11828. 



A. Targets 

A.l. r] Chamaeleontis Targets 

The i] Chamaeleontis cluster is locat ed at a distance of 97 pc. Age estimates of its 
constituent stars range fro m 4-10 Myr (ILawson et al.l 120011 ) . It is unusual compared to 
other star forming regions; ISicilia-Aguilar et al. J2009h measured a disk fraction of ~50%, 
much higher than expected for an older region. There is also a higher percen tage of disk 



systems in transition objects than found in other clusters ( IMegeath et al 



20051 ) 



A. 1.1. RECX 11 



RECX 11 is a weakly accreting CTTS. Its low accretion rate (M = 4xlO -11 M ?/r~ 1 ) 
would imply that ph otoevaporation should have cleared the inner dis k, but this is not see n 
( Gautier et al.l 120081) . UV-H2 emission implies the presence of gas ( IFrance et al.l 1201 lbl ). 
Ingleby et al.l (120111 ) also find low accretion values (M = 3xlO _lo M |/r _1 ) using the same 
data set presented here and argue that negligible planet formation has occured during the 
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disk evolution. From H2 line profiles, they measure rotational velocities fast enough to place 
the molecular g clS clS close as 3 R*. 



A.1.2. RECX 15 



RECX 15 is an active accre tor; the Ha EW is 90.0 A, and its profile is indicative of 
a blue-shifted wind component (ILawson et al.ll2004j ). The SED implies a continuous disk 

lar surface out to tens of AU or more, where the disk flares strongly 



from close to the ste 



(IGautier et a" 



2008 



Woitke et al.l l201lh . How ever, the H 2 (1 - 0) S(l) emission line width 



Ramsay Howat fc Greaved (120071 ) places the UV-H 2 a few AU from the star. 



measured by 

RECX 15 has a stellar mass of 0.2 M and inclination angle of 40° (jWoitke et al.l 1201 if ) . 



A. 2. Taurus- Auriga Targets 



The Taurus- Auriga complex is one of the nearest star forming regions at a distance of 
~140 pc. Its paucity of O/B stars, combined wit h a low interstellar radiation field, enable 
the persistence of many disks ( iLuhman et al.ll2010l ). Disk fractions in the Taurus region have 
been measured from ~ 45% for low-mass objects (M=0.01-0.3 M ) to ~75% for solar-mass 
stars. Some of these objects have been studied in exceptional detail, while others have only 
been measured in large surveys. 



A. 2.1. DM Tau 



DM Tau is considered to be a transitional disk due to its flux deficit shortward of 8 



/im ( jCalvet et al.ll2005l ). Resolved images and spectra of dust continuum emission show an 
inner hole of 19 AU with little dust inside this region (0.0007 lunar masse s) , although with a 
detectable accretion rate (~1.3xlO _8 M ?/r _1 ) there must be gas present (jCalvet et al.ll2005 



Andrews et al.l 120111 ) . The disk exhibits a rich chemistry; organic mol ecules CO, CN, CS , 
HCN, HNC, HCO + , C2H, and N 2 H + have been measured in emission (IDutrey et al .11199 71 ) . 
The abundances of these molecules show that stellar UV photons have strong influence over 
their c reation and destruction. DM Tau has a stell ar mass of 0.43 M and inclination angle 
of 35° jHartmann et al.lll998[ Andrews et al.lhoilh . 
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A. 2.2. HN Tau 



HN Tau is a binary system (0.7 and 0.4 M m) with strong jets that expand with distance 
from the primary star (IHartigan et al.l 120041 ) . In this work we study only the spectrum 
of HN Tau A, due to its separation from HN Tau B (3.1"). M-band observations from 



NIRSPEC on Keck reveal weak CO emission, but strong absorption 



communication). While its exact inclination is not known. France et al. 



J.M. B rown, private 
(J2011aj) measure CO 



electronic transitions detected in absorption, which could be indicative of disk gas viewed at 
high inclination. 



A. 2. 3. DR Tau 



DR Tau is a K4 star with an accretion rate that appears to vary by over two orders 
of magnitude ( iGullbring et al.l |2000| ; IJohns-Krull fc Gaffordl 120021 ). DR Tau shows strong 
molecular emission at IR wavelengths. Rotational transition emission lines of H oO are seen, 
in hi gh Js (20-50), as well as CO and OH at temperatures around 1000 K (jSalyk et al. 
20081 ). By modeling the gas as an isothermal LTE ring, the emission radii are 1.0, 0.8, and 
0.7 AU for H 2 0, CO, and OH, respectively. Time-lapsed CO meas urements show variability 
between overtone emission and absorption (ILorenzetti et al.ll2009l ). 



A. 2.4. DK Tau 



DK Tau is a binary system consisting of a K9 and M0 with ^ 350 AU separation . From 
mm imaging both appear to have individual circumstellar disks ( jJensen fc Akesonl 120031 ) , 
however our observation was targeted on DK Tau A. Based on the UV Unlike the majority 
of binary systems in which both stars contain disks, the A and B disks may be misaligned 
(IKraus fc Hillenbrand! 120091 ). Some molec ular emission has been detected in the DK Tau 
disks, particularly sub-mm 13 CO emission (jGreavesll2005l ). 



A.2.5. DF Tau 



A binary pair (M l + M3.5) with 13 AU separation, DF Tau may be as young as 0.1 
Myr (IGhez et al.lll997l ). Our spectrum samples the primary, as the secondary is very dim at 
UV wavelengths. Strong; a ccretion has been measured from a UV continuum ( IHerczeg et al. 
20061 ; France et al.l 1201 lbl ) . as have UV H 2 emission lines. H 2 appears both in Lya/C IV 
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pumped emission, and in absorption in an outflow ( lYang et al.l 1201 11 ). DF Tau displays a 
rich spectrum of IR-CO emission li nes; the CO (1 - 0) R(3) show s a FW HM of ~ 74 km s _1 , 



implying a close orbit to the star ( INajita et al. 



20031 ). iGreaved (120051 ) also detect sub-mm 



emission from 13 CO. DF Tau has a mass of 0.27 M and inclination of 85° ( iGullbring et al. 



1998l : IJohns-Krull & Valentil 120011 ) 



A. 2.6. LkCal5 



LkCal5 is a pre-transitional CT TS with a ringed gap in its disk, likely caused by a 
(proto)planet ( IKraus &: Ireland! 120111 ). From mid-IR imaging the inner rim of the optically 



thick dust disk is calcu l ated to be 46 AU, with some remnant dust inside (jEspaillat et al. 



2007; 


Pietu et al. 


2006 


)• 



0.1 AU. 



A. 2.7. GM Aur 



GM Aur is a massive transitional disk; the Spitzer Infrared Spectrograph (IRS) SED is 
best explained by 0.02 l unar masses of dus t within a 5 AU inner disk, and an outer disk with 
a 24 AU inner radius ( iCalvet et al.l 120051 ). High resolution sub-mm observations suppor t 
this picture , findin g a 28 AU hole in the optically thick dust disk ([Andrews et al.l 120111 ). 



Salyk et al.l (120071 ) measured near-IR CO emission within 0.21q2 AU, at high rotational 



temperatures, coincident with the measured dust inner rim. Ne II was m easured in the 
surrounding disk at the systemic velocity of 14.8 km s _1 (INajita et al.l 120091 ). Along with a 
profile broadened by rotational velocity, this would suggest the gas lies in the disk. 



A. 3. Isolated Targets 



Two of the targets do not reside in star-forming regions. V4046 Sagittarius is a ~9 Myr 
old CTTS, a closely orbi ting (~ 9 R separation ) binary with combine d mass of 1.78 M 
(jStempels k, Gahml 120041 ) at a distance of 83 pc (IRodriguez et al.l 120101) . At an inclinatio n 
of 36°, the disk app ears to be coplanar with the central binary stars (IMonin et al.l 120071 ). 
Kastner et al.l (120081 ) measured molecular mm emission lines from CO , HCN, CN, and HCO + ; 
the double peaked CO profiles confirm outer disk Keplerian rotation. IRodriguez et al.l (120101 ) 
imaged the system in 12 CO(2-l) and 13 CO(2-l) millimeter emission, measuring the molecular 
disk extent to ~ 370 AU and mass of ~ 110 Earth masses. HD 135344B is an F star in a 
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visual binary, and has a disk with a gap from .45 to 45 AU (IBrown et all 120071 ). It is more 
massive than the ot her stars in our sa mple. CO has been detected in emission inside the 
hole out to ~15 AU (IBrown et al.ll2009l ). The distribution o f gas column densities agr ee with 
models of a Jovian planet orbiting between 10 and 20 AU (IPontoppidan et all 120081 ) . 
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Fig. 1. — Extinction-corrected Lya profiles (left column) with geocoronal Lya removed, and 
CO (14 - 3) emission bands (right) for each of the sample targets. In the CO bands, an 
H 2 line is identified with a vertical dotted red line, and the dotted grey lines bound the 
high and low- J sides of the (14 - 3) band. The horizontal black/blue/red line identifies the 
local continuum, with the dashed orange lines at the la error in the continuum average. 
Representative error bars are shown in green. 
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Fig. 2. — Extinction-corrected Lya profiles (left column) with geocoronal Lya removed, and 
CO (14 - 3) emission bands (right) for each of the sample targets. In CO bands, an H2 line is 
identified with a vertical dotted red line, and the dotted grey lines bound the high and low- J 
sides of the (14 - 3) band. The horizontal black/blue/red line identifies the local continuum, 
with the dashed orange lines at the la error in the continuum average. Representative error 
bars are shown in green. DR Tau and LkCa 15 are non-detections. 
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Fig. 3. — (14-4) transitions for targets as identified. The horizontal lines show continuum 
levels, with dashed orange lines indicating the upper and lower limits for summation. The 
solid blue and red lines mark the wavelength range for the high and low- J regions as de- 
scribed in the text; vertical gray lines bound these regions. Vertical dashed red lines mark 
contaminating lines. Representative error bars are shown in green. H2 and CI I emission lines 
are identified with red-dotted lines, and O I emission features are identified in the spectra 
in which they appear. 
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Fig. 4. — V4046 Sgr Lya data with overlaid model profiles (top), and drawing of disk that 
produces the flux we see (bottom). Lya emission (red profile) is created near the star, then 
absorbed by an outflow (magenta profile), and additionally absorbed by interstellar H I 
(green profile) to produce the observed data. The CO in the disk sees the outflow-absorbed 
emission. Disk height is not to scale relative to extent. 
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Fig. 5. — DF Tau model results. Top: Lyct data overplotted with minimum, average, and 
maximum model Lya profiles. Middle: CO (14 - 3) band with model from average Tqo and 
Nco- Bottom: Tqo vs. Nqo for best fit models. The red square shows the values used in 
the middle plot. 
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Fig. 6. — Model-derived Tqo vs. literature accretion rates. 
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Fig. 7. — Ratios of blue to red flux for Lyct and CO. The blue and red Lyct fluxes were 
summed from 1210 - 1214.7 A and 1216.7 - 1222 A, respectively. 
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Fig. 8. — CO (14 - v") band luminosity vs. literature accretion rates for positive detection 
targets. The CO luminosities were calculated by applying branching ratios to the (14 - 3) 
luminosities. 
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Fig. 9. — CO (14 - v") band luminosities, divided by literature accretion rates, demonstrate 
an increasing correlation with age. 
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Fig. 10. — Top: RECX 11 (14 - 12) model emission (red) over the observed emisison. The 
data have been smoothed to one spectral resolution element for display purposes. Summed 
fluxes for the 1712 A emission vs. the 1316 A band, for each target. The dashed red line 
shows the expected (14 - 12) flux expected from branching ratios for a given (14 - 3) flux. 



Table 1: HST-COS observing log 



Object 




R.A. 


Dec. 


Date 


Texv (S) 


Program 


Object 




(J2000) 


(J2000) 




(G130M,G160M) 




RECX 11 


08 h 


47 m 


01.28 s 


-78° 59' 34.1" 


2009 December 12 


3645, 4514 


11616 


DF Tau 


0A h 27 m 


02.81 s 


+25° 42' 22.1" 


2010 January 11 


4828, 5549 


11533 


DK Tau 




44.25 s 


+26° 01' 24.5" 


2010 February 4 


1585, 1971 


11616 


LkCal5 


0A h 


39 m 


17.79 s 


+22° 21' 03.1" 


2010 February 5 


1712, 1725 


11616 


RECX 15 


08 h 


43 m 


18.43 s 


-79° 05' 17.7" 


2010 February 5 


3891, 4502 


11616 


DR Tau 


0A h 




06.22 s 


+16° 58' 42.6" 


2010 February 10 


1704, 1746 


11616 


HN Tau 


0A h 


33 m 


39.37 s 


+17° 51' 52.1" 


2010 February 10 


5725, 4529 


11616 


HD135344B 


15 h 


15 m 


48.42 s 


-37° 09' 16.3" 


2010 March 14 


2666, 4836 


11828 


V4046 Sgr 


18 h 


14 m 


10.49 s 


-32° 47' 34.2" 


2010 April 27 


4504, 5581 


11533 


GM Aur 


0A h 


55 m 


10.98 s 


+30° 21' 59.1" 


2010 August 19 


2128, 1861 


11616 


DE Tau 


0A h 


21 m 


55.69 s 


+27° 55' 06.1" 


2010 August 20 


2068, 1851 


11616 


DM Tau 


0A h 


33 m 


48.74 s 


+18 ° 10' 09.7" 


2010 August 22 


3459, 3770 


11616 



Table 2: Target Parameters 



Object 


bp. Type 


Group 


1 

d 


Ay 


M 


logio(Age) 


Reierenc 








(pc) 




(10 S M yr L ) 


(yrs) 




V4046 Sgr 


K5, K7 


Isolated 


83 





1.3 


7.08 ± 0.20 


6,8,10 


DF Tau 


Ml, M3.5 


Taurus-Aur 


140 


0.5 


17.69 


5.74 ± 0.15 


1, 11 


RECX 15 


M3.4 


rj Cha 


97 





0.1 


6.78 ± 0.08 


5,12 


DM Tau 


Ml 


Taurus-Aur 


140 


0.5 


1.12 


6.76 ± 0.20 


3,7,11 


RECX 11 


K6.5 


rj Cha 


97 





0.03 


6.60 ± 0.20 


5,13 


HD135344B 


F4 


Isolated 


140 


0.3 


0.54 


6.90 ± 0.30 


4,14 


DE Tau 


MO 


Taurus-Aur 


140 


0.62 


2.64 


5.82 ± 0.20 


1,11 


GM Aur 


K7 


Taurus-Aur 


140 


0.96 


0.3 


6.58 ± 0.25 


1,11 


HN Tau 


K5, M4 


Taurus-Aur 


140 


0.65 


0.13 


6.27 ± 0.24 


1,11 


DR Tau 


K4 


Taurus-Aur 


140 


1.2 


3.16 


5.29 ± 0.34 


2,9,11 


DK Tau 


K9, MO 


Taurus-Aur 


140 


1.42 


3.79 


6.17 ± 0.22 


1,11 


LkCal5 


K5 


Taurus-Aur 


140 


.62 


0.13 


6.35 ± 0.26 


3,7,11 



(1) iGullbring et all dl998h: (2) iGullbring et all (12 000) 



fl2006h: r5)lLawson et all (120041): (6) IStempels fc Gahm 



: (3) iKenyon fc Hartmannl (119951) : (4) iGarcia Lopez et al. 



2004h: m lHartmann et all (119981): f8l 



(120071 ): (13) lLawson et all (120011 ): (14) Ivan Boekel et all (120051 ) 



France et al. 



(201JJJ); 



(9) Iwhite fc Ghezl (boOlh : f lO) buast et al l j200ol): f Tlj jKraus fc Hillenbrandl J2009I): (12) iRamsav Howat fc Greaves 



-35 - 



Table 3: CO Fourth Positive bands identified in sample CTTSs. 



Band ID 


X obs (A) a Pumping Source 


Xpump 


(A) 


A(v' - v") (s- 


-i)b B{v'-v") c 


(14 - 2) 


1280.5 


H I Lya 


1214- 


1222 


4.0x10 s 


0.04 


(14 - 3) 


1315.7 


H I Lya 


1214 - 


1222 


1.5xl0 9 


0.17 


(14 - 4) 


1352.4 


H I Lya 


1214 - 


1222 


2.7xl0 9 


0.30 


(14 - 5) 


1390.7 


H I Lya 


1214 - 


1222 


1.6xl0 9 


0.17 


(14 - 7) 


1472.6 


H I Lya 


1214 - 


1222 


9.5x10 s 


0.10 


(14 - 8) 


1516.3 


H I Lya 


1214 - 


1222 


5.1x10 s 


0.06 


(14 - 10) 


1610.1 


H I Lya 


1214 - 


1222 


5.6x10 s 


0.06 


(14- 12) 


1713.2 


H I Lya 


1214 - 


1222 


3.9x10 s 


0.04 


(0-1) 


1597.3 


C IV 


1544 - 


1551 


4.3xl0 9 


0.33 


(0-2) 


1653.2 


C IV 


1544 - 


1551 


3.6xl0 9 


0.27 


(0-3) 


1712.4 


C IV 


1544 - 


1551 


1.8xl0 9 


0.14 



a CO wavelengths are taken from Kurucz (1993). 
''Summed A values for each (v' - v") band. 

c Branching ratios for each vibrational band, relative to all (v' - v") bands. 



Table 4: Best-fit Lya parameters 





li(0) 


°\ 


v 2 


i 2 (o) 


0"2 




Object 


(FEFU) 


(km s _1 ) 


(km s -1 ) 


(FEFU) 


(km s _1 ) 


(FEFU) 


V4046 Sgr 


10122 ± 3531 


263 ± 40 


-72 ± 35 


5232 ± 2840 


368 ± 31 


7.0 ± 3.0 


DF Tau 


52084 ± 943 


374 ± 








4.5 


RECX 15 


16 ±» 


26^ + 292 

zuo -265 


208 ± 28 


906 ± 146 


224 ± 10 


0.8 ± 0.2 


DM Tau 


is±fi 


339+11 


-133 ± 20 


270 ± 64 


437 ± 14 


4.0 


RECX 11 


1376 ± 1 


197 ± 








2.0 


HD 135344B 


2487 ± 1109 


346 ± 52 


118 ± 151 


2293 ± 2168 


158 ± 95 


0.8 ± 0.2 


DE Tau 


4835 ± 4590 


249 ± 59 


-13 ± 78 


lfiOQ+1914 
iuyo -1693 


225 ± 152 


1.1 ± 0.5 



Columns: (2) Peak intensity of rest gaussian emission component; (3) Width parameter for rest gaussian emission 
component; (4) Peak intensity of secondary gaussian emission component; (5) Width parameter for secondary gaussian 
emission component; (6) Local continuum for Lya 



Table 5: Best-fit Lya parameters (continued) 





Vout 


N out 






Object 


(km s^ 1 ) 


(10 19 cm- 2 ) 


(10 19 cm^ 2 ) 




V4046 Sgr 


-96 ± 15 


3.0 ± 0.8 


1.1 ± 0.4 


104 - 136 


DF Tau 


-36 ± 7 


14.0 ± 1.5 


49.0 ± 1.4 


< 1.08 


RECX 15 


-85 ± 36 


1.1 ± 1.0 


0.9 ± 0.2 


2.2 - 20 


DM Tau 


119 ± 178 


1.4 ± 1.2 


41.8 ± 2.1 


< 1.16 


RECX 11 


-99 ± 4 


1.9 ± 0.1 


3.3 ± 0.1 


< 13 


HD 135344B 


-99 ± 26 


5.6 ± 2.7 


7.1 ± 1.6 


0.76 - 7.1 


DE Tau 


-103 ± 56 


11.4 ± 7.7 


34.8 ± 8.3 


< 1.11 



Columns: (2) Outflow absorber velocity; (3) Outflow absorber column density; (4) ISM absorber column density 



Table 6: Best-fit CO parameters 









Continuum 


^^broad 


log(N C o) 


Tco 


vco 




Object 


Band 


(FEFU) 


(km s -1 ) 


(cm -2 ) 


(K) 


(km s -1 ) 






V4046 Sgr 


14 - 


3 


3.1 


27 ± 1 


18.55 ± 0.45 


241 ± 41 


-2.6 


2.3£ 


\ - 5.2 


V4046 Sgr 


14 - 


4 


3.5 


30 ± 3 


18.47 ± 0.41 


295 ± 87 


-6.4 


5.18 


- 7.55 


DF Tau 


14 - 


3 


2.8 


21 ± 1 


18.56 ± .38 


434 ± 135 


27 


< 


1.08 


RECX 15 


14- 


3 


0.4 


16 ± 1 


18.50 ± 0.53 


318 ± 131 


1.0 


< 


1.14 


RECX 15 


14 - 


4 


0.5 


20 ± 1 


18.54 ± 0.57 


427 ± 249 


1.0 


1.7 


- 2.04 


DM Tau 


14 - 


3 


0.5 


19 ± 4 


19.21 ± 1.14 


833 ± 562 


35.0 


< 


1.18 


DM Tau 


14 - 


4 


0.8 


17 ± 3 


19.08 ± 0.77 


775 ± 445 


26.2 


< 


1.14 


RECX 11 


14 - 


3 


0.2 


53 ± 33 


18.36 ± 0.50 


347 ± 158 


18.0 


< 


1.13 


HD 135344B 


14 - 


3 


1.0 


41 ± 14 


18.58 ± 0.93 


148 ± 66 


5.0 


< 


1.13 


DE Tau 


14 - 


3 


0.3 


26 ± 20 


18.35 ± 1.03 


793 ± 696 


30.0 


< 


1.14 



Columns: (3) Local continuum for CO band; (4) FWHM of individual J emission lines; (5) CO column density; (6) 
CO temperature; (7) CO LSR velocity 
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Table 7: CO Band Summed Fluxes 



Target 


C(14-3) 


Blue F(i4_ 3 ) 


Red F(i4- 3 ) 


L(CO) 




(FEFU) 


(FEFU A) a 


(FEFU A) 


(10 28 ergs s- 1 ) 


V4046 


3.30 ± 0.68 


47.36 ± 0.39 


20.33 ± 0.28 


33.70 ± 0.24 


DFTAU 


3.73 ± 1.02 


4.69 ± 0.31 


17.38 ± 0.40 


31.25 ± 0.71 


RECX15 


0.54 ± 0.31 


2.11 ± 0.12 


3.35 ± 0.15 


3.71 ± 0.13 


DMTAU 


2.20 ± 1.02 


9.19 ± 0.41 


4.21 ± 0.39 


18.98 ± 0.80 


RECX11 


0.33 ± 0.24 


2.18 ± 0.11 


1.01 ± 0.10 


2.17 ± 0.10 


HD135344B 


1.28 ± 0.73 


6.38 ± 0.39 


3.18 ± 0.35 


13.55 ± 0.75 


DETAU 


1.54 ± 2.24 


3.04 ± 0.52 


0.99 ± 0.49 


5.70 ± 1.01 


GMAUR 


2.82 ± 1.37 


2.70 ± 0.40 


0.60 ± 0.39 


4.68 ± 0.79 


HNTAU 


2.16 ± 1.12 


1.30 ± 0.29 


1.43 ± 0.31 


3.87 ± 0.60 


DRTAU 


14.79 ± 8.81 











DKTAU 


4.14 ± 4.90 


4.57 ± 1.49 


8.16 ± 1.74 


18.02 ± 3.25 


LkCal5 


1.61 ± 1.50 












a l FEFU A = 1CT 15 ergs cm- 2 s" 1 



